CAFPE-139/10 
FTUV/10-0428 
IFIC/10-10 
UG-FT-269/10 



O 
(N 



in 



00 



Pinched weights and DuaUty Violation in QCD Sum Rules: 

a critical analysis 



^ ■ Martin Gonzalez- Alonso'^, Antonio Pich" and Joaquim Prades^ 



^ ! " Departament de Fisica Teorica and IFIC, Universitat de Valencia-CSIC, 

> : Apt. Correus 22085, E-46071 Valencia, Spain 

CAFPE and Departamento de Ffsica Teorica y del Cosmos, 



G^ ■ Universidad de Granada, Campus de Fuente Nueva, E- 18002 Granada, Spain 

^' 

o 
o 



X 
}^ 

d • Abstract 

We analyze the so-called pinched weights, that are generally thought to reduce 
the violation of quark-hadron duality in Finite-Energy Sum Rules. After showing 
how this is not true in general, we explain how to address this question for the 
LR correlator and any particular pinched weight taking advantage of our previous 
work [Tj, where the possible high-energy behavior of the LR spectral function was 
studied. In particular we show that the use of pinched weights allows to determine 
with high accuracy the dimension six and eight contributions in the operator product 
expansion, Oq = (-4.3^°;?) • lO'^ GeV^ and Og = i-7.2tti) ■ GeV^ 



1 Introduction 



In a recent work [Ij we have analyzed the violation of quark-hadron duality (DV) of a given 

t(0+1) / 
'■ud,LR\ 



QCD Sum Rule with the non-strange LR correlator n(g^) = ^udLRi^"^) defined by 



KlLuis) = ^y"dSe^'^^(0|T(LC,(x)K,(0)t)|0) 



= i-g^'q' + q'qn ^^fniq') + g^'^q' nS,^(g^) , (1) 

where L^^i^) = m7''(1 - 75)6? and R^di^) = u'y^{l + 75)^. 

A QCD Sum Rule [2] takes advantage of the analytic properties of the correlator to 
relate its imaginary part in the positive real g^-axis (where hadrons lie) with its value in 
the rest of the complex plane, where the Operator-Product Expansion (OPE) allows us 
to calculate it in terms of quarks and gluons: n*-'™(s) = Xlfc ^2fc/(— s)'^. The DV comes 
from the fact that this OPE breaks down in the vicinity of the positive real g^-axis. We 
can write a general QCD Sum Rule for the LR correlator in the following form 

[ " ds w{s) p{s) + — (f dsw{s)U^^^{s) + DV[w(s),so] 

= 2/Xm2) + Res[«;(s)n(s)], (2) 

where p(s) = ^ Imn(s) and w(s) is an arbitrary weight function that is analytic in the 
whole complex plane except in the origin (where it can have poles). The violation of 
quark-hadron duality is formally defined as [T1|3H7] 

DV[^(.),.o] = ^ / d. ti;(.) (n(.) - n«p^(.)) . (3) 



Using analyticity one can write the DV in the following form 

/•oo 

I)V[w{s),so] = / dsw{s) p{s) , (4) 



so 



that shows how the DV is nothing but the part of the integral of the spectral function 
that we are not including in the sum rule. In Ref. we have studied the DV from this 
perspective, using the following parametrization 

p{s>Sz) = K e~^^ sm{f3{s - Sz)) , (5) 

for the spectral function beyond ^ 2.1 GeV^ and finding the region in the 4-dimensional 
parameter space that is compatible with the most recent experimental data [9J and the 
following theoretical constraints: first and second Weinberg Sum Rules [10] (WSRs) and 
the sum rule of Das et al. [TT] that gives the electromagnetic mass difference of pions (vrSR). 
The parametrization emerges naturally in a resonance-based model [21[I2l[I3] that has 
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been used recently to study the violation of quark-hadron duality [SJ [TUTU [15], although 
without imposing the previously explained theoretical constraints in the numerical analysis. 

In Ref. [T] we used this parametrization to calculate the DV associated to Finite- Energy 
Sum Rules (FESRs) with the weights w{s) = {n = —2, —1, +2, +3), but it can be used 
to analyze any other QCD Sum Rule with the LR correlator. In this letter we would like to 
apply the results of [I] to the so-called pinched-weight FESRs, where the standard weight 
s" is substituted by a polynomial weight that vanishes at s = Sq (or near this point). 

It has been often assumed that the use of pinched-weights (PWs) minimizes the DV0 
[^l5l [T6]423] ■ since they suppress the contribution from the most problematic region in the 
contour integral of Eq. ([3]) , close to the real axis [21] . However the alternative expression for 
the DV given in Eq. @ shows that things are more subtle [HEllT] and that the assumption 
is not necessarily true, since a PW function will indeed suppress the first part of this 
hadronic integral but at the same time may enhance the high-energy tail that can become 
important. If the final balance is positive and the weight function does its job minimizing 
the DV contribution is something that depends on the particular weight used and on how 
fast the spectral function goes to zero, something that is not known theoretically. 

This question about the convenience of the use of these PWs is very entangled with the 
more general question of how to estimate the duality violation of a given sum rule. The 
observation of a more stable plateau in the final part of the data range is the standard 
requirement to check if the weight improves the situation, and the deviations from the 
plateau the standard way of estimating the remaining DV. However it is important to 
notice that the existence of the plateau is a necessary but not sufficient condition, because 
it could be temporary. This is particularly plausible because the PWs produce curves that 
have derivative zero in the second duality poinlH (sq ~ 2.6 GeV^), which is very near of 
the end of the data. That is, they produce a fake plateau, that can induce to the possibly 
wrong conclusion that the DV is negligible for that weight and that value of Sq. 

Here fake means that the correlations between the experimental points of the plateau 
are extremely high and such that we do not have several points indicating the same value, 
but just one point drawn several times. In principle, a fit of these points to a straight 
line is sensitive to the correlations and would tell us if the plateau is real or it has been 
artificially created by the weight function, but in practice this is not always possible, since 
the high correlations among points prevent us from using the standard x^"fit, as explained 
in [250. 

The results obtained in our previous analysis [Ij allow to address these issues in a 
quantitative way. In particular, we will study the PW versions of the QCD Sum Rules of 
Ref. pp, i.e. those that give an estimation of the hadronic parameters Cg^, L^q, and Og. 

^It must be emphasized that the PW functions are also useful because they are expected to minimize 
the experimental errors, since they suppress the region near the kinematical end point. 

^Thc duality points are two particular points (located at sq 1.5 GeV^ and sq 2.6 GeV^) where 
both the first and second WSRs happen to be satisfied, i.e. where their DV contributions vanish. 

•^This situation was found e.g. in the determination of the V-A condensates in Ref. [4l[5]. 
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2 Numerical analysis 



We are interested in PW functions that do not introduce new unknown quantities (con- 
densates of higher dimension), since in that case a clean analysis is not possible anymore, 
and more specifically we will work with pinched- weights w{s) that have a double zero in 



s — ^pwi thSiij IS 



Sth V ''pw / \ fpw J Opw ^pw 



ds^(l-—Y = -8Lf, - 4^ + 2^ - DV[t._i, sj , (7) 



Sth V ^pw / Spw 



[ ^ds p{s){s-sp^f = 2fls%-Aflmlsp^ + 2/>t + O, - DV[w2, s.] , (8) 

>^Sth 



Sth 



ds p(s) (s - Spwf {s + 2s 



pw ) 



Sth 



6/>^s' + 4/^^^ + 2/>^ - - DV[^3, . (9) 



The results depend on the point Spu, where the weight is pinched. In order to suppress 
the experimental error it is convenient to pinch the weight at the left of the matching point 
Sx, whereas in order to suppress the DV-error (dispersion of the histograms) it is convenient 
to pinch it at the right of s^. We have scanned the region finding that the optimal choice 
of Spw, that is, where the errors are minimizecill, is Spw ^ Sz ^ 2.1 GeV^. 

A careful comparison between these PWs and the standard weights s"" shows that in 
the case of the condensates the former are smaller (in absolute value) than the later for 
any s > Sp^, and therefore are expected to generate a smaller DV0, although the question 
of how large is the remaining DV is not clear at all. In the case of the chiral parameters 
L^Q and Cgy, the convenience of the PW is not known a priori and it depends essentially 
on how fast the spectral function goes to zero, in order to suppress the enhancement that 
the PWs produce in the high-energy region (see Eq. (jlj)). In other words, the key point is 
the value of the 7 parameter, that is around one [1]. We will show that this value is large 
enough to suppress the high-energy tail and so to benefit from the use of the PW. 

In Ref. yy we have used the parametrization ([5]) for the spectral function p(s) and 
we have analyzed the allowed parameter space once the experimental and theoretical con- 
straints are taken into account. In other words, we have generated a large number of 
"acceptable" spectral functions, compatible with both QCD and the data. The differences 
among them determine how much freedom is left for the behavior of the spectral function 
beyond the kinematical end of the r data. In particular we can calculate the value of the 
parameters Cg^, L^q, Oq and Cg obtained through the sum rules ([S]) - ([2]) for each of these 



Obviously the optimal point is different for every sum rule ([6]) - ([9]), but the difFerences are negligible 
within errors. 

^Notice that this is not a mathematical statement, but only a hand- waving estimate and can be altered 
due to accidental cancellations. 
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Figure 1: Statistical distribution of values of C^^ (upper-left), L'jq (upper-right), Oq (lower- 
left) and Os (lower-right) for the accepted spectral functions, using the pinched-weight sum 
rules (ED - ([9D with Spin — ~ 2.1 GeV . The parameters are expressed in GeV to the 
corresponding power. 



possible spectral functions. The results of this process are given in Fig. [H which shows the 
statistical distribution of the generated values. We can see that the histograms are much 
more peaked around their central values than those obtained in Ref. [1] with standard 
weights. 

Let us remind that in addition to the error associated to the DV (estimated from the 
dispersion of the histograms) we have the experimental ALEPH error, and both depend 
on the used weight function. In principle one expects the PWs to minimize also the 
experimental uncertainties, since they suppress the region near the kinematical end poin10. 

The associated numerical results are (we give the 68% probability region). 

(8.168^[!:S[}f ±0.12) -10"^ GeV~2 = (8.17 ± 0.12) -10"^ GeV-^ (10) 

(-6.444 ± 0.05) ■ 10"^ = (-6.44 ± 0.05) ■ 10"^ , (11) 

{-4.33 toil ± 0.65) ■ 10-^ GeV^ = (-4.31°;?) ■ 10^^ GeV*' , (12) 

(-7.2 tli ± 2.9) • 10-=^ GeV^ = (-7.2 l^J) ■ 10"=^ GeV* , (13) 

where the first error is that associated to the high-energy region (integral from to 
infinity) , that we compute from the dispersion of the histograms of Fig. [H and the second 
error is that associated to the low-energy region (integral from zero to Sz), that we compute 

^Notice below that in the case of Os this does not happen. This is because the PW enhances the 
low-energy region errors sizably. 



/-left 

T eff 
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in a standard way from the ALEPH data. For the sake of comparison we show the analogous 
results obtained in Ref. [1]: Cg = (8.17 ± 0.12) • lO'^ GeV'^ Lfl = (-6.46 t [J:[]?) ■ 10-^ 
Oe = {-5Atli) ■ 10-3 GeV^ and = {-8.9t\Y) ' lO'^ GeV^ where we can clearly see 
the improvement achieved with the PWs. 

Since the first error in Eqs. (fT0l) - (fT3l) is not Gaussian we show also here the 95% prob- 
ability results: 

C^^ = (8.168^°;°°^ ±0.24) -10-3 GeV-2 = (8.17 ±0.24) -10-=^ GeV-^ (14) 
= (-6.444 ±0.1) • 10^3 = (-6.4 ±0.1) ■ 10-3 , (15) 
Ce = (-4.33+^:™ ± 1.3) ■ 10"^ GeV^ = (-4.31?;^) • 10"^ GeV*' , (16) 
Os = [-7.21:^,^ ± 5.8) ■ 10-3 GeV^ = (-7.2 1^2^) " IQ-^ GeV^ . (17) 



3 Beyond the dimension eight condensate 

We can play the same game with higher- dimensional condensates, where using again 
pinched weights w{s) that have a double zero in s = Spw we have 



/ ds p(s) (s - Sp^f (s^ ± 2sp^s ± 3s5^) 



-8/>2s3 + Qfy + 2/>8 + 0,0 - DVK, s.] , (18) 



/ ds p(s) (s - Sp^f (s3 + 2sp^s^ ± 3Sp^s ± 4Sp^) 



-10/>^,4 ^ 8^2^5 ^ 2/>f - Ou - DVK, s.] , (19) 



/ ds p(s) (s - Spujf (s^ ± 2sp^s3 + 3Sp^s^ ± 4s3^s ± 5Sp^) 



/" ds p(s) (s - Sp^)' (s^ ± 2sp^s^ ± 3s;^s3 + 4s3^s2 ± 5sJ^s ± 6sJ^) 

= -14/>2,6^ + i2/>J^ ± 2/>i^ - - DV[wr, s.] . (21) 

Working again with Sp^ ~ s^ ~ 2.1 GeV^ we find the results shown in Fig. |2l The 
associated numerical values are (68% C.L.) 

O^o = (±4.11}:^)- 10-2 GeV^V (22) 

Oi2 = (-0.12l°:°3^) GeV'\ (23) 

= (±0.2 «i) GeV'\ (24) 

Ci6 = (-0.2 l°i) GeV^^ , (25) 

where all the errors come from the dispersion of our histograms since the experimental 
error is very much smaller for these higher- dimensional condensates. The 95% probability 



5 



tt tuples tt tuples 




Figure 2: Statistical distribution of values of Oio, 12,14,16 for the accepted spectral functions, 
using the PW sum rules f|T8l) - fl2T|) with Spw = Sz ~ 2.1 GeV^. The parameters are 
expressed in GeV to the corresponding power. 



results are: 

Oio = (+4.1 tlD • 10"^ GeV^° , (26) 

O12 = (-0.121°:}^) GeV^^ , (27) 

Ci4 = (+0.2 ± 0.5) GeV^^ , (28) 

C16 = (-0.2 1};?) GeV'^ . (29) 

It is really impressive that the sign of the condensates can be established for Oiq and O12 
since the importance of the high-energy region in their determination is huge. One could 
have expected that the differences between our possible spectral functions would generate a 
huge error in these higher- dimensional condensates, but our conditions (WSRs+vrSR+data) 
have turned out to be very restrictive about the acceptable spectral functions allowing quite 
precise extractions. 



4 Comparisons and summary 

We have used the method developed in Ref. [Ij to analyze the error of different pinched- 
weight Finite-Energy Sum Rules and to extract the value of different hadronic parameters. 
Comparing the results obtained here with those of Ref. [I] we see that, as theoretically 
expected, the use of the pinched weights is less beneficial to the determination of the low- 
energy constants L^q and Cg^ than to the determination of the condensates. Our final 
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Figure 3: Comparison of our results for Oq (left) and Os (right) with previous determina- 
tions [H[5|[^[T ^ [2U |l7rHiD] (we show for every method the most recent determination). The 
blue bands show our results at 65% C.L., while the 95% probability regions are indicated 
by the dotted lines. 



results for the former are in excellent agreement with the most precise determination of 
them [26]: Cgf = (8.18 ± 0.14) ■ 10^^ GeV"^ and Lf^ = -(6.48 ± 0.06) • 10"^ Notice 
that, even if these determinations are also based on the PW sum rules ([6]) and ([7]), the 
estimation of the error presented here is obtained through a completely different method, 
based on more solid grounds and represents a confirmation of them. 

We have obtained quite precise measurements for the condensates Oq and Os using the 
PW sum rules ([H]) and ([9]). In this way we have checked that the PW succeeds in minimizing 
the errors and we can conclude that the most recent experimental data provided by ALEPH, 
together with the theoretical constraints (WSRs and vrSR), fix with accuracy the value of 

and almost determin^ the sign of Oq. Our results are compared in Fig. [3] with previous 
determinations of Oq and Oq. One recognizes in the figure the existence of two groups of 
results that disagree between them. For Oq there is a small tension between a bigger or 
smaller value, whereas in the case of Og the disagreement affects the sign and is more 
sizable. As can be seen in Table [H these discrepancies also appear in higher- dimensional 
condensates, that we have also extracted applying the same method. 

Our results agree with those of Ref. [HElIiniEO] since they also use pinched weights, 
but we think ours are based on much more solid grounds, due to the completely different 
approach followed. We see in fact that the DV error was underestimated in Refs. [ll|5], 
especially in the determination of the higher-dimensional condensate^. 

^One can see in our final result (fT7|) that at 2a a positive value of Os is already allowed, but it must not 
be forgotten that the distribution is highly non-Gaussian and we can see in the corresponding histogram 
of Fig. [1] that the possibility of being positive is negligible. 

^This is just an explicit case where we can see that even when the pinched weights generate less DV 
than the standard weights s", the observed plateau is in part artificially created and hides the DV. That 
is why the errors of Ref. [5] are underestimated. 
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Oio X 10^ 


Oi2 X 10^ 


Cl4 X 10^ 


OiQ X 10^ 


This work 


' —Id 


-120+^!] 


±200+^Jl[i 


-200+^[!!] 

— 4UU 


A/Tpi cii n PIT! Rj Ppri rTQl 

IVJ-CXO J LLCLll Oo J- v>l lO lO 1 


—14 ± 12 








Narison [22] 


-17.1 ±4.4 


±14.7 ±3.7 


-9.6 ±3.1 


±4.3 ± 1.9 


Friot et al. [37J 


-13.2 ±3.6 


±13.3 ±3.9 


-12.8 ±3.9 


±11.9 ±3.8 


Zyablyuk [33j 


-4.5 ±3.4 








Almasy et al. |10] 


±66^^!^ 








Bordes et al. [20] 


±72 ± 28 


-240 ± 50 






Latorre & Rojo [M] 


±78 ± 24 


-260 ± 80 






Cirigliano et al. [5] 


±48 ± 10 


-160 ±30 


±430 ± 60 


-1030 ± 140 



Table 1: Comparison of our determination of 0io,i2,i4.i6 with other works. The condensates 
are expressed in GeV to the corresponding power. The results shown for Ref. [5] are 
those obtained with the old ALEPH data (with the OPAL data the numbers are not very 
different), and the results shown for Ref. [37] are those obtained with the minimal hadronic 
ansatz, that is, without the addition of the p' resonance, that in any case modifies just 
slightly the results. 

We also agree with the results of Ref. [30] based on the use of the second duality point, 
although that technique has a much larger error. It is remarkable also the agreement with 
Ref. 03] that is the only one that follows a technique similar to ours, trying to analyze the 
possible behavior of the spectral function but through a neural-network approach. Their 
result has a bigger uncertainty, maybe only due to the fact they used the old ALEPH data. 

Our analysis indicates that the DV error associated to the use of the first duality 
point is very large and was grossly underestimated in Ref. [331 [35], where also higher- 
dimensional condensates were neglected. In Ref. [37 t [39 | l^ the numerical values obtained 
at this first duality point are supported through theoretical analyses based on the so-called 
"minimal hadronic ansatz" (a large- A'"(7-inspired 3-pole model) or Pade approximants. Our 
results show however that the first duality point is very unstable when we change from the 
WSRs to the C6,...,i6 sum rules, indicating that the systematic error of these approaches 
is non- negligible. Essentially the same can be said about Refs. [28||29] where the last 
available point sq = rn^ was used. The minimal hadronic ansatz [37l|41] gives a reasonable 
approximation to O^, but its accuracy seems not good enough to reproduce the signs of 
the higher-order condensates (although an alternating-sign series is indeed predicted). 

Summarizing, our results agree within two sigmas with the other estimates of Oq, but 
for condensates of higher dimension 08,io,i2,i4,i6 they agree with Refs. [Hl5|[T9 | 120 1 130 1 131] 
but not with Refs. [28 ] |29 | |33 | [351437] . It is worth noting that in particular our method 
shows that and are both negative, whereas it suggests that the sign alternates for 
higher-dimensional condensates. 
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